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| Bianchi type IX, 'Mixmaster' universes are investigated in low-energy- 



effective-action string cosmology. We show that, unlike in general relativity, 
there is no chaos in these string cosmologies for the case of the tree-level 
action. The characteristic Mixmaster evolution through a series of Kasner 
epochs is studied in detail. In the Einstein frame an infinite sequence of 
chaotic oscillations of the scale factors on approach to the initial singularity 

X" 

S-h ' is impossible, as it was in general relativistic Mixmaster universes in the 

presence of a massless scalar field. A finite sequence of oscillations of the scale 
factors described by approximate Kasner metrics is possible, but it always 
ceases when all expansion rates become positive. In the string frame the 
evolution through Kasner epochs changes to a new form which reflects the 
duality symmetry of the theory. Again, we show that chaotic oscillations 
must end after a finite time. The need for duality symmetry appears to 
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be incompatible with the presence of chaotic behaviour as t — > 0. We also 
obtain our results using the Hamiltonian qualitative cosmological picture for 
Mixmaster models. We also prove that a time- independent pseudoscalar axion 
field h is not admitted by the Bianchi IX geometry. 
PACS number(s): 98.80.Hw, 04.50.+h, 11.25.Mj, 98.80. Cq 
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I. INTRODUCTION 



String cosmology has attracted a lot of interest recently, especially in the context of 
duality symmetry, which is a striking feature of the underlying string theory. The motion 
of the simplest bosonic string in background fields is governed by the nonlinear sigma- 
model action Cosmological solutions (to lowest order in a' - the inverse string tension) 
have been considered in detail in many papers with special interest in the possibility of 
inflation, the behaviour of inhomogeneities, and the relation between the so called pre- and 
post-Big-Bang phases of evolution. M. The bosonic string spectrum of particles contains 
the graviton, dilaton and axion (antisymmetric tensor field). The dilaton can always be 
accommodated within homogeneous geometries, but this is not the case for the axion. It 
has been shown that the antisymmetric tensor-field potential does not necessarily have 
to be time-independent and this can prevent the isotropisation of homogeneous models 
at late times ||. Moreover, the time- dependent antisymmetric tensor potential cannot 
be accommodated in all homogeneous anisotropic universes of Bianchi or Kantowski-Sachs 
types |3J]. In this paper we concentrate on Bianchi IX (BIX) axion-dilaton cosmology with 
homogeneous ansatz as given originally in |J . This is equivalent to the inclusion of a time- 
dependent pseudoscalar axion field, h = h(t). The dilaton field is always homogeneous, 
= 0(t). 

It is well known that in the vacuum BIX homogeneous cosmology one approaches the 
initial singularity chaotically ||. An infinite number of oscillations of the orthogonal scale 
factors occurs in general on any finite interval of proper time including the singularity at 
t = 0. These oscillations are created by the 3-curvature anisotropy of the spacetime and 
are intrinsically general relativistic in origin. Physically, the propagation of homogeneous 
gravitational waves alters the curvature of spacetime along the direction of propagation, so 
that their non-linear back-reaction on the curvature reverses the direction of propagation. 

In the presence of a massless scalar field (or, alternatively, 'stiff matter', with pressure, p, 
equal to density, p) the situation changes. Only a finite number of spacetime oscillations can 
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occur before the evolution is changed into a state in which all directions shrink monotonically 
to zero as the curvature singularity is reached and the oscillatory behaviour ceases |7[. In 
this paper we want to investigate the oscillatory approach to singularity in BIX low-energy- 
effective-action homogeneous string cosmologies with both axion and dilaton fields. 

This investigation complements other recent studies of the behaviour of string cosmolo- 
gies at early times. The investigation in || shows that there is an open neighbourhood of 
initial data space for string cosmologies which display velocity-dominated behaviour. That 
is, on approach to the singularity the spatial gradients become negligible compared to time 
derivatives, velocities are non-relativistic, and curvature anisotropies are negligible. These 
solutions resemble inhomogeneously-varying Kasner metrics, but the approximation only 
hold on scales larger than the particle horizon. Barrow and Kunze || found a large family 
of exact inhomogeneous string cosmologies which are velocity-dominated. They describe the 
evolution of the dilaton and axion fields on all scales, showing how these fields oscillate once 
inhomogeneities in their distribution enter the horizon. These studies show that the general 
behaviour of the low-energy-string cosmology equations are considerably simpler than those 
of general relativity. The string gravity sector behaves in a quasi-Newtonian fashion in the 
vicinity of cosmological singularities. In general relativity the velocity-dominated, quasi- 
Kasner, behaviour is disrupted by the chaotic oscillations of the three-curvature anisotropy 
on approach to the singularity whenever the matter content is a massless scalar field or 
perfect fluid satisfying p < p. It is therefore important to determine what happens in the 
presence of three-curvature anisotropies in string cosmology. It is possible that there exists 
another stable neighbourhood of evolution on approach to the initial singularity which is 
chaotic rather than of quasi-Kasner form. We will investigate this situation by studying the 
Bianchi type IX string cosmologies in order to discover whether chaotic behaviour is possible 
or whether it is destroyed by the presence of the dilaton and axion fields. 

The paper is organized as follows. In Section II we derive the string low-energy-effective- 
action equations for BIX geometry. We present the equations in both the string and the 
Einstein frames and we use orthonormal frames in the main body of the paper. In Section 
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Ill we discuss the possible routes to chaos near a curvature singularity in both the Einstein 
and the string frames with reference to the calculations given in Ref. ||. In Section IV we 
discuss the relation between the oscillating Mixmaster behaviour of the scale factors and the 
duality symmetry of the string equations. In Section V we use the Hamiltonian approach 
to represent the evolution of the BIX string cosmology as the motion of a 'universe point' 
inside a curvature potential. This picture provides a simple picture of the conditions needed 
for chaos to occur, and for it to be absent, in the Einstein, string and so-called axion frames. 
In Appendix A we give the relation between Ricci tensors for BIX model in orthonormal and 
holonomic (coordinate) frames and show why it is impossible to include a time-independent 
pseudoscalar axion field h in the non-axisymmetric and axisymmetric BIX models. 

II. LOW-ENERGY-EFFECTIVE-ACTION EQUATIONS FOR BIANCHI IX 

UNIVERSES 

The low-energy-effective-action field equations in the string frame are given by |3],[5] 1 

R v » + v M V0 - \n mP E^ = o, (2.i) 

R - V^W + 2V M V> - ^H^pH^ = 0, (2.2) 

V M ( e -*iP" a ) = dp (e^^H^) = 0, (2.3) 

where <fi is the dilaton field, H^p = Qd\p,B v p\ is the field strength of the antisymmetric tensor 
B^ v = —B Ufl . In the Einstein frame we have || 

fy-\%R = J(?p+ ) +fp*)), (2.4) 

1 Following standard notation we use Greek indices here to write down the field equations (2.1)- 
(2.8). However, beginning with the formula (2.9) we replace these Greek indices by Latin ones 
i,j, k, I = 0, 1, 2, 3 because we are using an orthonormal basis. Greek indices a, (3,fj,,v = 0, 1, 2, 3 
are coordinate basis indices and will be used in Appendix A 
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V M V^+V 2< ^ 2 = 0, (2.5) 
6 

V M (e~ 2<t> H^ a ) = ~d, (e-^yf^H^ = 0, (2.6) 

and (k 2 = 8ttG, c = 1) 

K 2fuW = 1 f ~A~„ CT _ 1 ~,~a.\ v A 0V a 0, (2.7) 



2 2 JM ' . 

k 2 t;^ = im m ,H^ - k# 2 ) • (2.8) 
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Covariant derivatives are formed with respect to the Bianchi type IX metric which, in the 
string frame, reads as [[TJJ 

where the orthonormal forms cr 1 ,^ 2 ,^ 3 are given by (i,j,k,l = 0,1,2,3 are orthonormal 
basis indices) 

a 1 = cos ipd9 + sin ip sin 9 dip, (2-10) 

o" 2 = sin ipd9 — cos ip sin 9dp, (2- 11) 

a 3 = dtp + cos 9d<p, (2.12) 

and the angular coordinates ip, 9, <p span the following ranges, 

< ^ < 4tt, < 9 < 7r, < ip < 2tt. (2.13) 

In the Einstein frame the metric is 

rf5 2 = d p _ ~ a \t){a 1 ) 2 _ b 2 (t)(a 2 ) 2 - ~c 2 (t)(a 3 ) 2 , (2-14) 

and 

dt = e'^dt, (2.15) 

Oi = e-^ 2 ai, (2.16) 

where a, = {a, 6, c, } and a; = {a, 6, c}. 



The nonzero Ricci tensor components in the orthonormal frame a 1 are (an overdot means 
a derivative with respect to the synchronous coordinate time t) B 



- Rq 



- R 3 
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a b 




a c 
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+ 
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a b 




a c 
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a b 




b c 




+ 




+ 




" b 


a b 


b c 


c 




a c 




be 




+ 




+ 




c 




a c 




b c 



b 2 -c 2 



(a 2 -c 2 )' 
(a 2 -b 2 )' 



and the Ricci scalar reads 
-R 



a b c ab at be 
2- + 2- + 2- + 2-r + 2— + 2 7 - 

a b c ab ac be 

-TT^T-M + b 4 + c 4 - 2a 2 b 2 - 2b 2 c 2 



2a 2 c 2 



(2.17) 
(2.18) 
(2.19) 
(2.20) 



(2.21) 



Since the model under consideration is homogeneous the dilaton field can only depend on 
time and so we have (cf. Appendix A) 



V V°^ 
V o 0V°< 



V) 

di 



(2.22) 
(2.23) 
(2.24) 



As for the axion field, we can use the following homogeneous ansatz |||| 
H = —e 1 A e 2 A e 3 = Aa 1 A a 2 A a 3 , 



abc 



where A is constant and 



(2.25) 



e = a,;(j . 



[2.26) 



Since 



(Hi) 2 = H ikl W k \ 



(2.27) 



and 
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Tj2 TT zjikl 



(2.28) 



we have 



0\2 



my = {m 



H 2 



2\2 



[H. 



3\2 



-2 



A 2 



a 2 b 2 c 2 



A 2 



a 2 b 2 c 2 



(2.29) 
(2.30) 

(2.31) 



The special case of axial symmetry a = b gives the orthonormal forms (2.10)-(2.12) as 
below fll[ 



a 1 = dO, 
a 2 = sin 

a 3 = dip + cos 



(2.32) 
(2.33) 
(2.34) 



One can easily check (cf. Appendix A) that the homogeneous ansatz (2.25) is also appro- 
priate for the axisymmetric case. 

III. TIME-DEPENDENT PSEUDOSCALAR AXION FIELD AND THE 
OSCILLATORY APPROACH TO A SINGULARITY 



The ansatz (2.25) is, in fact, equivalent to the inclusion of a time-dependent pseudoscalar 
axion field h = h(t) || (cf. Appendix A). We are going to investigate the possible emergence 
of chaos in such models. Using (2.17)-(2.24) and (2.29)-(2.31) the field equations (2.1) in 
the string frame read (' = d/dt) 



d db dc d ■ 

- + -- + (J> 

a ao a c a 

b db be b ■ 

b ab be b 



2a 2 b 2 c 2 



2a 2 b 2 c 2 
1 



— a 



c ^ dc _|_ c c ; 
c ac be c 2a 2 b 2 c 2 



(b 2 -c 2 ) 
(a 2 - 
! Ir) 



d b c 

- + T + -- 

a b c 

1 A 2 



c z ) -b q 

2 



2 a?b 2 c 2 
1 A 2 
2a 2 b 2 c 2 
1 A 2 
2a 2 b 2 c 2 



0. 
0. 
0. 
0. 



(3.1) 
(3.2) 
(3.3) 
(3.4) 



S 



Eq. (2.2), with Eq. (3.1), reads 



2a 2 b 2 c 2 



fab ac bc\ -fa b c\ - 9 
\ab ac be J \a b c) 

1 V + 6 4 + c 4 - 2a 2 6 2 - 2aV - 26V 1 ' ^ 



+ 



2 a 2 b 2 c 2 



0. 



(3.5) 



Using (3.5), together with the sum of (3.2)-(3.4), we have the dilaton equation of motion 
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b" 




c" 
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a" 




a' V 


a' 
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a 




a b 
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b" 




a! V 
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c' 




+ 








b 




a b 


c 


c 


c" 




c'b' 


a' 


~d 




+ 








c 




c b 


a 


c 



V 2 

T 

1 



!/~2 



i A2 


= 0. 




(3.6) 


a 2 b 2 c 2 




become 


(' = 


d/dt) 




A 2 \ , 
a 2 b 2 c 2 






(3.7) 


\b 2 -~c 2 ) 2 


-a 4 " 




(3.8) 


"(a 2 -c 2 ) 2 


-b 4 ' 




(3.9) 


\a 2 -b 2 ) 2 


-c 4 " 




(3.10) 



and (2.5) now becomes 



\a b c 



-2(b 



A 2 



a 2 b 2 c 2 

A new time coordinate is introduced to simplify the field equations by defining 

1 ~ 

drj = ——dt = —=r-dt. 

abc dbc 

First, notice that the string-frame equation (3.6) then simplifies to 

(j) m - A 2 e- 2<t> = 0, 



(3.11) 



(3.12) 



(3.13) 



where (, v = d/drj). The Einstein-frame Eq.(3.11), with the time coordinate t, gives the same 
result, (3.13). The solution of (3.13) is || 



= cosh AMr] + \ 1 



A 2 



A 2 M 2 



sinh A Mr], 



(3.14) 



with M, A constant (A 2 M 2 > A 2 ). A useful relation, implied by (3.14) is 
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0, w + 0^ = A 2 M 2 



(3.15) 



Let us introduce new forms for the scale factors 



a = e 



b = & 



c = e ' 



and 



a = e 



b = e? 



c = e ' 



so, from (2.16) 



a = a — 0/2 



/3 = /3-0/2 



7 = 7- 0/2, 



(3.16) 



(3.17) 



(3.18) 



and dt = exp (<j)/2)dt in (3.12). The field equations (3.1)-(3.4) in the string frame take the 
form 



(a + (3 + 7) jJW — M 2 = 2 («„/?, + a ,7, + /3,,7,»j) ~ 2 («,»? + A»j + l, v ) <f>/ v , 
2e 2<t> a vv = (b 2 -c 2 ) 2 -a 4 + A 2 , 
2e 2 ^, w = (a 2 -c 2 ) 2 -6 4 + A 2 , 
2e 2 ^7, w = (a 2 -6 2 ) 2 -c 4 + A 2 . 

The equations (3.20)-(3.22), using (3.13), can be rewritten to give 



(-<!>+ m, m 

("0 + 27) irm 



(b 2 -c 2 ) 2 -a 4 
(a 2 -c 2 ) 2 -b 4 
(a 2 -b 2 ) 2 -c 4 



-20 



-20 



-20 



(3.19) 
(3.20) 
(3.21) 
(3.22) 



(3.23) 
(3.24) 
(3.25) 



Notice that there is no explicit dependence of the axion, A, in these equations (3.23)-(3.25). 
On the other hand, using (3.12), Eqs.(3.7)-(3.10) in the Einstein frame reduce to 



2 + a v % + p„<y tV ) = (a + (3 + 7) + ^M 2 , 



>w 2 



2a 



2/3 )TW = (a 2 - g 2 



(ga-c 2 ) -a 4 , 

27, w = (5 2 -6 2 ) 2 -c 4 . 
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(3.26) 
(3.27) 
(3.28) 
(3.29) 



These equations (except for the constant M which appears in (3.26)) have exactly the same 
form as the standard BIX equations of general relativity ||. Equations (3.26)-(3.29) can be 
explicitly transformed into (3.23)-(3.25) by using (3.16)-(3.18). 

In order to discuss possible emergence of chaos we consider suitable initial conditions 
expressed in terms of the Kasner parameters. 



A. Einstein Frame 

The Kasner solutions are obtained as approximate solutions of the equations (3. 7)- (3. 11) 
when the right-hand sides (describing the curvature anisotropics) are neglected. In the 
Einstein frame, in terms of i-time, they are 



a = a t Pl , 



b = b F\ 



c = c W\ 



while (A 2 < A 2 M 2 ) 



for A = 0, and 



In d + In 



for A 0, where 



~ M ,~. 
In a + — In \t , 
A 



A 2 \ ,~, M 

= t~T 



AM 2 



(3.30) 



(3.31) 



(3.32) 



A = a b c Q , 



d = const. 



(3.33) 



From (3.8)-(3.10), irrespective to the presence of the axion term A, we have the following 
algebraic conditions on the Kasner indices, pi : 



(3.34) 
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and, from (3.7), 



Jt, 9 M 2 . , 

= 1 3.35 

f={ 2A 2 



which is exactly the case considered first by Belinskii and Khalatnikov [J/J So, in the 
Einstein frame their argument follows. 
In terms of 7/-time from (3.27)-(3.29) 

a = Kpit] + f\, 

(3 = Ap 2 r] + f a , (3.36) 
7 = Ap 3 r] + f 3 , 

and the constraint equation (3.26) bocomes f\ 

1M 2 . . 

PXP2 + PXP3 + P2P3 = ~j ■ (3.37) 

This constraint is equivalent to the constraints (3.34)-(3.35). In order to get (3.30) from 
(3.36) we need to use (3.12) to relate t-time with ^-time, i.e., 

r] = A" 1 In t + const. (3.38) 

When we approach singularity i] — > — oo (t — > 0) we cannot neglect all the terms on the 
right-hand side of (3.27)-(3.29) since one of the Kasner indices in (3.37) can be negative. 
We assume p\ = —\px\ < 0, which means that a{rf) ^> b{rf) ^> c{rj), so (3.27)-(3.29) are 
approximated by 



2 In their notation it happens if we put M 2 /2A 2 = q 2 . Also, they define the Lagrangian of the 
scalar field without a factor of 1/2 (i.e., = (V</>) 2 ) while we follow standard notation. 

3 In || there is ^4 2 instead of M 2 while in all our calculations we obtained M 2 both in the string 
frame and in the Einstein frame. They claim they use the so-called 'sigma frame', which, however, 
seems to be just the Einstein frame [fl0|. 
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1 4a 

2 




1 4d 


l>m = 


2 



id 



(3.39) 



Far away from singularity the approximate axisymmetric Kasner regime is fulfilled, but 
it is broken by the a 2 term when we approach singularity, so our Kasner solutions (3.36) will 
be fulfilled for 77 — > 00 (t — > 00); suitable solutions of (3.39) which fulfil (3.36) are 



5(77) = 


-iln 

2 




Iln" 


Hv) = 


2 


7(77) = 


Iln' 
2 



— - cosh (— 2 1 pi I A77 
2|pi|A V IWI 

— - cosh f— 2 1 »i I A77 
2|pi|A V mi ' 



+ A(-|pi| +p 2 )f/, 
+ A(-|pi| +p3)v- 



In the limit 77 — > —00 (t — > 0) from (3.40) we have 



(3.40) 



5(77) 


= A|pi|77 


= -Api77, 




/%) 


= A (pa - 


2|pi|)r/ = 


A(p 2 + 2^1)77, 


T(»7) 


= A (pa - 


2\pi\)v = 


A (p 3 + 2p\) 7), 



(3.41) 



which means that one Kasner epoch (3.36) with indices pi is replaced by another Kasner 
epoch with indices given by (3.41). By virtue of (3.12) 



77 = A' hit + const., 



(3.42) 



and 



4 This necessarily requires p~\ < 0. Of course, one has the same solutions for p\ > 0, but in such a 
case our initial conditions would have to be taken at rj — > 00 rather than at r/ — > —00. This does 
not change the physics of the problem and was assumed in [19|,21|, but we prefer to follow the spirit 
of 0. 
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b = b' t p z, (3.43) 



where 



N 


Pi 


Pl ~ l-2|px| 


l + 2pa : 


-/ P2-2|pi| 


_ p 2 + 2pi 


P2 " 1-2N 


1 + 2pi ' 


-/ Pa-2|pi| 


P3 + 2pi 


% " l + 2fr 


1 + 2pi ' 


A't = abed, 


A' 

X ~ 



(3.44) 



and 



3 3 1 M 2 

Ep5 = L Epf = l" ^- (3-45) 

i=i i=i ^ A' 2 

As we can see from (3.32), the axion does not influence these asymptotic solutions. 
One can easily show that 

[2 M [2 , t \ 

~ V3 - X7i " V3 < 3 ' 46 > 

and, after ordering the Kasner indices by pi < p2 < ps, we require 

- 1/3 < Pi < 1/3, 

< p 2 < 2/3, (3.47) 
l/3<p 3 <l. 

This means that, unlike the vacuum case, all the Kasner indices can be positive and, as in the 
analysis of 0, the final situation is that the universe inevitably reaches a monotonic stage of 
evolution in which all three Kasner indices are positive (see also the discussion of Section V)Q. 



5 Then, the terms of the type e 4a ,e 4 ^,e 4 ^ decrease if r\ — > —00 (t — * 0), and they do not allow 
a transition into another Kasner epoch to take place. This is just monotonic evolution like in 
isotropic Friedman case. 
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By means of Eqs. (3.30) and (3.32), which lead to the same relations between the Kasner 
indices (3.34)-(3.35), one can extend this conclusion into the axion-dilaton cosmology. The 
presence of the axion field cannot change the fate of the universe near to a singularity and the 
universe finally reaches a monotonic stage of evolution with all three scale factors tending to 
zero as t — > 0. Thus, there is no chaos in BIX homogeneous string cosmology in the Einstein 
frame. Finally, we note that there is an isotropic limit here, provided all the Kasner indices 
are equal with pi = p 2 = p3 = 1/3. 



B. String Frame 



For the string frame we can follow the discussion given in [|T^] (see also [ 20|,2T1 ). 



In terms of t-time, Kasner solutions of the system (3.1)-(3.4), for A = 0, are given by 

a = a t Pl , 

b = b t p \ (3.48) 
c = c t P3 , 
e "t> = d t p \ 

or, alternatively, the last of these conditions (3.48) reads as 

<f>(t) = -In d -P4ln|t|, and P4 = -M. (3.49) 
After putting (3.48) into (3.6) we have 

p 4 = l-f>, (3.50) 

i=i 

and from (3.1) we have 

i> 2 =i, ( 3 - 51 ) 

1=1 

so (3.49) can be rewritten to give 

3 

0(t) = -lndo + (5>i-l)ht- (3-52) 



i=i 
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The last condition means, for instance, that a choice p, = (—1/3,2/3,2/3) gives J2Pi — 1 
and leads to a vanishing or constant dilaton (M = 0), i.e., to general relativity. (Other 
permutations of the three pi with the signs changed are allowed). Thus, the difference 
between general relativity and string theory depends on the "fourth Kasner index" p^ = —M. 
Now, from (3.23)-(3.25) (one can see that Api = 2kp i here) 

a = 0/2 + A(pi/2)77 + ri/2 = 0/2 + Ap^ + f u 

(3 = 0/2 + A(p 2 /2)r] + r 2 /2 = 0/2 + Ap 2 r] + f 2 , (3.53) 
7 = 0/2 + A(p 3 /2)r] + r 3 /2 = 0/2 + Ap 3 r] + f 3 , 

and from (3.19) 

M 2 

+ pm + P2P3 = -j2~- (3.54) 

The condition (3.54), in fact, can be obtained from (3.51)-(3.52). In order to get (3.48) 
we need to use (3.12), i.e., 

r] = A -1 hit + const., A = a boCod . (3.55) 

From the conditions (3.50)-(3.51) on the Pi, we have necessarily that 

-l-V3<M<-l + V3. (3.56) 

However, the domain of M given by (3.57) covers the whole duality-related region since in 
the case of a Kasner regime the duality symmetry (cf. Section IV) simply means that we 
change 

pi -> -pi ,p 2 -^-p2 ,P3-*~P3 ,M -> M -pi -p 2 -p 3 = -{M + 2). (3.57) 

Having given (3.58), we can delineate the two duality-related domains of Kasner indices as 
follows: 

_ 1 _V3<M<-l,--L< Pl < v | ,_!<„,< ,_ 1 < P3 <_I 1 | ; (3.58) 
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and 

-iSW<-i + ^,-^|<p.<^ ,-\^<v,<\ \<^<v,<i; (3.59) 

Of course, for —1 — < M < — 1, the indices are ordered so that p^ < P2 < Pi; while, for 
— 1 < M < — 1 + v^3, they are ordered as p\ < p 2 < p%- 

From (3.57)-(3.59), we can draw some interesting conclusions. First, that the vacuum 
general relativity case M = (with —1/3 < p% < 0,0 < P2 < 2/3,2/3 < p% < 1 ) is dual 
to the case M = -2 (with -1 < p 1 < -2/3, -2/3 < p 2 < 0, < p 3 < 1/3), while the case 
M = — 1 is 'self-dual' in M giving M — — 1 again, although the p'^s change. Second, the 
M = — 1 plane is the dividing plane for the duality- related range of the parameters which 
are defined by 

l/2 1/2 
-2V3 - P2 - 2V3' 

The isotropic Friedmann solution is given when we choose M = — 1 + y/3, p\ = P2 = 
P3 = while its dual with M = —1 — y/3, pi = p 2 = ps = —1/V3 does not appear 

in the vacuum general relativity case. This shows that one can generally have two different 
types of change in the Kasner indices, and in M (effectively, a fourth Kasner index ). There 
are oscillations as in general relativistic vacuum or stiff fluid cases [f||7]], but there are also 
duality-related exchanges of indices. We will discuss some points related to duality that 
refer to the work of ref. |T6| . 

There are eight possible permutations of the Kasner indices in which some are equal 
to —l/y/3 or 1/a/3, with suitable M. These are the first two quadruples (pi,p 2 ,p 3 ,M) = 
(1/V3, l/v 7 ^, l/y/3, -1 + v/3), (-l/v 7 ^, — 1/V3, — I/a/3, -1 - >/3), phis another three pairs 
with (-1/V3, 1/V3, 1/V3, -1 + (V3/3)), -l/v 7 ^, -l/v 7 ^, -1 - (V3/3)) and the per- 

mutations which are dual to each other. One can easily show that these last three pairs 
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of cases are duality related and describe exactly the transitions from one Kasner epoch 
to another given by (3.80). Another interesting set of 'self-dual' (M = —1) combina- 
tions which describe LRS (Locally Rotationally Symmetric) Kasner solutions are given 



by (^2/3, -1/y/E, -1/V6, -1), (- V 2 A l/VQ, -1) and their permutations. Other 

interesting duality-related combinations are (—1/3,2/3,2/3,0), (1/3,-2/3,-2/3,-2) and 
(0, l/y/2, -l/y/2, -1), (0, -l/y/2, 1/V2, -1). The last two give the so-called Taub points 
(flat Minkowski spacetime) ||16|| . 

The question arises whether the changes of the Kasner indices (chaotic and duality- 
related) have something in common. In order to answer this question one should try to 
find a suitable parametrization of the Kasner indices analogous to the one given in J7j for 
the stiff-fluid case. We have checked that such a parametrization does not cover the right 
range of the indices given by (3.58) unless we multiply p\ by y/3. After some searching, we 
have found the following parametrization of the indices which has advantages which we will 
explain in due course. 

The w-parametrization we use is defined as follows f\ 

Pi = T~r~~~2 ' ( 3 - 6 °) 

1 + u z 



1 1 

Vi = - 



P 3 = ol~< — 2 
21 + u 

where u is constant and 



(l + M)(l + w 2 ) -2u + P + (M,u) 
(1 + M)(l + u 2 ) -2u- P + (M, u) 



21 + u 2 

11 . , 



P + (M, u) = V(l - 2M - M 2 ) (1 + it 2 ) + 4m (M + 1) (1 + u 2 ) - 12m 2 , (3.61) 
and to achieve pi < we need u < 0. Alternatively, we have 



6 In the general relativity case (M = 0) this parametrization can be written down in terms of the 
one suggested in |l6| with p\ = (1/3)(1 — 2 cos V0,£>2,3 = (1/3)(1 + cos-0 ± sin^) with a suitable 
choice of cos?/' = (1/2) (it 2 — 6u + 1)(^ 2 + 1) in our case. 
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ft 1 i 9 ' 

1 + U 



V2 = ~ 
P3 



' 1 '1 + M)(1 + u 2 ) +2u + P-{M, U ) 



21 + u 2 

11 .. 



21 + m 2 



(l + M)(l + u 2 ) + 2u- P~(M,u) 



where 

P-(M, U ) = \] {1 — 2M — M 2 ) (1 + u 2 ) 2 -4m (M+ 1) (1 + w 2 ) - 12m 2 , (3.63) 

and for p\ < we need u > 0. These are very general transformations that include all 
duality-related cases (3.58)-(3.59) and the vacuum general relativity case when M = 0. 
Similarly, as for stiff-fluid models in Ref. 0, we draw a plot the allowed values of the 
parameter u against M and show the duality-related regions in Figs.l and 2. 

There are many possible transformations of the general homographic type u — > (au + 
b)/(cu + d) or some more general types related to Pade or other rational approximants 
|T8| which may be useful in the discussion of Mixmaster oscillations. First, notice that a 
transformation of the type 

u — > — , (3.64) 
u 

gives 

Pi(M, l/u) = pi(M, u), (3.65) 
Pi (M,l/u)= Pi (M,u), 

where i = 1,2,3, unless we take negative value of the root of P ± after applying the trans- 
formation. On the other hand using 

u -»■ — , (3.66) 
u 

changes parametrization (3.60) into (3.62), i.e., 

Pi(M, -l/u) = Pi {M, u) = -Pi(M, u), (3.67) 
Pi (M, -l/u) = Pi{M, u) = - Pi (M, u). (3.68) 
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Second, one can consider the transformation 

u — > —u, (3.69) 

which results in the same rules as (3.65). However, this transformation is not enough to 
describe duality relations between the considered solutions. As one can see, the duality- 
reflecting transformations should be of the type (see also Fig. 1) 

u->-u, M->-(M + 2), (3.70) 

or, alternatively 

M->-(M + 2). (3.71) 

u 

From (3.70) and (3.71) we see that 

P ± (-M - 2, -u) = P ± (-M - 2, -l/u) = P ± (M, u), (3.72) 

and 

Pi(-w) = -Pi(u), (3-73) 
P2(-M-2,-u) = -p 3 (M,u), 
p 3 (-M-2,-u) = -p 2 {M,u); 

similarly, 

p 1 (-u) = -p 1 (u), (3.74) 
P2(-M-2,-u) = -p 3 (M,u), 
p 3 (-M-2,-u) = -p2(M,u). 

We notice that the duality in (3.73)-(3.74) would entail the exchange of the Kasner indices 
p 2 and p 3 . 

If we assume that a ^> b ^> c in (3.23)-(3.24) then we obtain the following set of 
approximating equations (for A = 0) 
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4oc-2AMr) 




P,rm = 


1 

2 


(3.75) 


7,w = 


1 4a -1AMr\ 

2 




(P,rm = 


0. 





together with the constraint (3.19). The Kasner solutions in terms of 7/-time are given by 

ot{rj) = Apii] + const, 
(3{rj) = Ap 2 r] + const, 

7(77) = Ap 3 r] + const, (3.76) 
(f)(r)) = AMt] + const. 

However, these solutions are not directly obtained by using the relations (3.18) between 
the Einstein frame and string-frame scale factors. This is very important in the case when 
the axion field is taken into account and will be discussed below. 

The solutions of equations (3.75) which fulfil the above conditions, (3.76), in the limit 
77 — > 00 (t — > 00 - i.e. far from singularity), provided pi — — \pi\ < and 2p x — M = 
—2\pi\ — M < 0, can be chosen to 



«0?) = ( A(2pi 1 _ M ) coshA ( 2 Pi - M )^j + ^ AM *7. 

0(V) = ( A ( 2pi 1 - M) C ° Sh A(2P1 ~ M)r/ ) + (Pl + - M) Ar; + ^AMr/, (3.77) 

7(77) = i In f ^^Z M) cosh A ( 2 Pi - M )^) + (Pi + Ps ~ M) A V + l -AM V . 

In the limit 77 — > —00 (t — > - i.e. on the approach of singularity) they approach the 
following asymptotic forms Q 



7 Note we can derive them in a similar way to the for A 7^ case (see the rest of this section) , but 



here it is more convenient to follow the results given in [21]. 
8 If we assume pi < and 2pi — M < 0, then p\ — M < 0, provided M > 0, which means we have 
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ot{ri) ~ -A(pi - M)?7, 

P(ri) ~ A(p 2 + 2pi - M)t7, (3.78) 
7(77) ~ A(p 3 + 2pi - M)t7, 
<p(rj) ~ AM?]. 

One can check the solutions (3.77) by putting them into the constraint (3.19) in order to 
recover the condition (3.54), as expected. 

Now, we can express the scale factors in terms of the new Kasner parameters 

a = a' t p>1 , 

b = b' t p i, (3.79) 
c = c' t?K 
e -+ = d' t p '*, 

where 

Pi-M 

Pi = 



and 



P2 



l + 2pi — M' 
p 2 + 2gi - M 
1 + 2pi - M ' 



= P3 + 2 Pl -M 

1 + 2p! — M v ; 

-.1/ 

P4= l + 2 Pl -M = - M - 



A' = a' b' c' d' , 

v = (A')- 1 In t + const., (3.81) 



A' = (1 + 2pi - M)A. 



changed expansion of a(r/) into contraction. 
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If we take the axion field into account (A ^ 0) and assume that a ^> b 3> c in (3.23)- 
(3.25)), then we obtainP] 



OL m = ~ I A- - t ) < 

P,m = ~ {A 2 + e 4 °) e~ 2 \ 

l,rm = 



- (A 2 + e 4a ) e~ 2 \ (3.82) 



2 

A 2 e- 2 * 



with 0(r/) given by (3.14). The Kasner solutions in terms of rj-time are now given by (compare 
(3.18)) 

a(rj) = - (pi + M) 1] + const = Agi?7 + const., 
P(rj) = - (j>2 + M) rj + const = Aq^r] + const., 

7(77) = — (ps + M) 7] + const = Aqii] + const., (3.83) 
4>(rj) = AMr] + const. 

The solutions which fulfil the above initial conditions (3.83) for rj —>■ oo (p± < 0) are 

a (v) = ~\ ln coshpx^ + ^0(77), (3.84) 

l/l \ 1 1 

(3{rj) = - ln I — cosher/ J + - (pi + p 2 ) 77 + -0(r?), (3.85) 

7(77) = i ln coshp!^ + ^ (p! + p 3 ) r] + ^(77), (3.86) 



0(?7) = ln 



/ A 2 
cosh Mr/ + A / 1 — sinh M77 



(3.87) 



9 If M > 0, then the term e -2fli increases for r\ — ► —00. If, in turn, pi < 0, M > 0, and 2pi — M > 0, 
then the term e 2 ( 2 Pi- M )v decreases for r\ — > —00 and the whole picture is dominated by the axion 
term l/2A 2 e~ 2 ^ . It follows that the field equations become isotropic a m = fi m = j lrrn = \j2$ m = 
l/2A 2 e~ 2 ^ . The conclusion is that axion isotropises the model and chaos is impossible in such a 
case. It seems that the time-dependent axion field ansatz (cf. Appendix A) would allow chaos, but 
it is not admitted by the BIX geometry. 
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or, alternatively 



yu 2 \2q 1 -M 



cosh(2 gi -M)77 +-0(77), 



1 



«")=2 h 2 91 -M 



1 \ 1 

cosh(2 ?1 -M)?7j + ( gi + g 2 - M) 77 + -^(77), 



7 ^ = 



cosh(2g 1 -M)r/j + ( 9l + g 3 - M) 77 + -^(77), 



(3.88) 
(3.89) 
(3.90) 



with 0(77) unchanged. 

One can easily check by putting these solutions into the constraint (3.19), that the condi- 
tion (3.54) is fulfilled, which in turn ensures that the conditions (3.50)-(3.51) are fulfilled. In 
particular, note that, for (3.88)-(3.90), we need to replace p^s by q[s. In the limit 77 — > —00, 
(that is t — > 0), they approach the following forms 

®(v) ^(Pi - M K 

/3(77) ~ ^(p 2 + 2pi-M)77, 
T(»7) ~ ^(p 3 + 2pi-M)?7, 
0(7/) ~ --M77, 



(3.91) 



or 



a(77) ~ -A(gi - M)r], 
(3(7]) rsj A(q 2 + 2q x - M)r], 
7(77) ~ A(g 3 + 2 Pl - M)r), 
0(77) ~ -AM77. 



(3.92) 



Having given the conditions (3.50)-(3.51), one can express the indices p 2 and p 3 by using p 1 
and M, i.e., 



P3= 2 



(M + 1- Pl ) 



-3pj + 2 Pl (M + 1) + 1 - M (M + 2) 



(M + 1 - Pl ) + yj-3pl + 2p 1 (M + 1) + 1 — M (M + 2) 



(3.93) 
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Since the expression under the square root should be nonnegative, one can extract the 
restriction (3.56) on the permissible values of M. However, we are interested in knowing 
whether the curvature terms on the right-hand side of the field equations (3.23)-(3.25) really 
increase as i] — > — oo (t — > 0). This would require either a 4 e~ 2< ^, b 4 e~ 2< ^, or c 4 e~ 2< ^ to increase 
if the transition to another Kasner epoch is to occur [ W ,^0]] ■ Since 

frV 2 ^ oc t {2p2 ~ M) = £( 1+ ^3-p0 (3.94) 

we need one of the following three conditions to be fulfilled (remember that we have assumed 
Pi < 0,M > 0): 

2pi - M = 1 + pi - p 2 - p 3 < 0, 

2p 2 - M = l+p 2 -p 3 - Pl <0, (3.95) 
2p 3 - M = 1 + p 3 - pi - p 2 < 0. 

The three conditions (3.95), with the help of (3.93), are equivalent to 

M 

Pi < y, (3.96) 
-3p 2 + 2pi (M + 1) + 1 — M (M + 2) > 0. (3.97) 

A plot of these conditions is given in Fig. 3. The last of these conditions, (3.97), provides 
bounds on the possible values of M if a transition to occur: 

-2<M<jj. (3.98) 

Now, we see that the regions where the Friedmann isotropic limit is possible (all the Kasner 
indices equal - this happens for M = —1 — and M = — 1 + y/3) are excluded. One can 
always find the range of the indices for a transition from one Kasner epoch to another to 
occur in the string frame. 
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Instead of expressing the conditions for Kasner-type transitions in terms of pi and M, 



we can follow the pattern of and write them in terms of pi and p 2 - From the conditions 
(3.50)-(3.51), we can write 

p 3 = M+l- Pl -p 2 , (3.99) 



M = Pl +p 2 -l± ^l-pl-pl (3.100) 

where the plus sign is for M > — l,p 3 > and minus sign for M < — l,p 3 < 0. So, pi and 
Pi must be such that (M real) 

\-p\-p\>^ (3.101) 
and one of the three conditions (3.95) must be fulfilled, i.e., either 

Pi + p\ - P1P2 + Pi - P2 < 0, 

P 2 i+pI-PiP2-Pi+P2<0, (3.102) 



for p 3 > 0, or 



p\ + vl + P1P2 ~ Pi - P2 > 0, 



Pi + P\ ~ PlP2 + Pi ~ P2 < 0, 



Pi +P2 -P1P2 ~Pi +P2 < 0, (3.103) 

Pl+P2+ PlP2 + Pi + P2 > 0, 

for p 3 < 0. 

The plot of these conditions is given in Fig. 4. 

In summary, we have been able to determine the range of values that can be taken by 
the Kasner indices in the string frame and have proposed a parametrization which describes 
the evolution of these indices. Also, we have determined the values of the Kasner indices 
(see Figs. 3 and 4) for which the spacetime oscillations can really take place. However, 
now we have to determine whether the oscillations can be stopped once they have started 
as t — > 0. Before we come to this in Section V we first discuss some relations between the 
Kasner indices and duality. 
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IV. EXCHANGE OF KASNER INDICES AND DUALITY 



The class of homogeneous solutions of the low-energy-effective-action equations of string 
cosmology have a global 0(3,3) invariance also called T-duality (other duality invariance 
called ^-duality also appears in superstring models, see e.g. [Q) under which (0 is the 
so-called shifted dilaton field) 

M — > M' = f2 T Mf2, = 0- In VdetG -> 0, (4.1) 

where O is 6 x 6 constant matrix satisfying 

o T no = n, n = (J J) > (42) 

and 1 is the 3x3 identity matrix and 

m 4g G -g"- G bg B 4 (4 ' 3) 

where G = and B = By are 3x3 matrices. Any 6x6 constant matrix O obeying 
Eq.(5.2) generates new solutions M' from the original set M. If the antisymmetric tensor 
field (axion) vanishes under a choice Q = II, we deal with the "scale factor duality" where 



a' 2 



a' 

b' 2 - ^ (4-4) 



1 

2' 



0' — > — 2 In a&c. 

It is useful to define the logarithm of an average scale factor f3 and the so-called shifted 
dilaton defined by fl5fl 

A = -/=lnaj, (4.5) 

/3 = -^=ln (abc), 
= 0- y^. 
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Using (3.40)-(3.45) we have the relations (4.5) in terms of Kasner indices, i.e., 

4>= 0-lnA- (M + l)lnt, 

fa = -j= (liiaoi +pdnt) , (4.6) 

P= -p [lnA + (M + l)lnt] , 
V3 

where a 0i = {a , b , c }. In these variables the duality symmetry is just expressed by 

ft(t) = -A(t), (4.7) 

P(t) = -P(t), 

fa) = fa); 

or, in terms of Kasner indices, it reads pi — > —pi, M + 1 —>■ — (M + 1). After inclusion of 



time symmetry we have [15 



A(t) = -Pi(-t), (4.8) 
fa) = -fat), 

In the isotropic case pi = ±l/v3 and we recover exactly the case given in Ref. [15 



The same relations can be written down using the time coordinate r\ instead of t. Using 
the exact expressions, (3.76), for Kasner solutions we have 

fi = A( M + 1)^ 

= - y/3p - a v , 

A = AfiV- (4-9) 
So, one could relate these by duality symmetries 

fa) -> fa V ), (4.10) 
P(V) - -P(-V), (4-11) 
for r] — > ±oo respectively, using chaotic changes pi — > — p*. 
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V. HAMILTON! AN APPROACH TO BIX STRING MODELS 



In this Section we formulate a generalized Kasner model in Hamiltonian formalism as in 



Refs. ||i7||19|p2| in order to discuss the conditions for an infinite sequence of scatterings to 
occur against the walls of the curvature potential. As in the previous sections we discuss 
the problem in both the Einstein and the string frames. We also introduce the so-called 
B-frame, or axion frame pi, in which axion is minimally coupled. 



A. Einstein Frame 



We introduce the following standard parametrization for the Einstein-frame scale factors 



a = e 



c = e 



(5.1) 



and we define the potential, which describes the spatial curvature anisotropy (2.21) felt by 
scale factors in the Einstein-frame by 



V(ip±) = e~ 2a V^ ± ), 



(5.2) 



where 



1 r 



e -8V>+ + 2e 4 ^+ (cosh 4^3^- - l) - 4e~ 2 ^+ cosh 2^- 



Using (5.1)-(5.2), the Eqs. (3.7)-(3.10) read as 



(5.3) 



a 



12 ~ ^ + V- + — 0' 2 + — A 2 e~ 2<t> - && + -e- 2& V(ij±). 



12 



12 



6 



(5.4) 



This is the Hamiltonian constraint. The Einstein-frame action in terms of the scale factors 
(5.1) after integrating out spatial variables is given by (compare f23| ) 



S = dte 



^3a 



-6a' 2 + 6^' 2 + 6C 2 + \<P' 2 + ^W 2 + e~ 2& V^ ± ) 



(5.5) 
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and the conjugate momenta are 



7T+ = 12^e 35 , 
tt_ = 12^j'_e 3& , 
vr = 0'e 33 , 

7r CT = a'e 2 ^ 35 = const. = A, (5.6) 

so the Hamiltonian is 

# = _!i + zi + !^ + !i + !i + e A& V(?P ± ). (5.7) 
24 24 24 2 2 K ^ ±J v ; 

Now, we follow the standard discussion of the potential walls 

V = e 45 V(^±) (5.8) 



being hit by a particle moving in the potential well (see Ref. [0]). In the region ip + <C — 1 
and ■?/>_ ~ 0, the approximate distance from the origin of coordinates ip + and to the wall 
is given by 

D = ~&, (5.9) 
while the maximum apparent velocity of this wall is 

Vmax = O l . (5.10) 

The velocity of a particle moving against the walls is 

v p = vV? + ^- 2 , (5-11) 

and it will not be scattered infinitely many times if there is some region of the potential 
which the particle enters and from which it cannot catch up with the wall, i.e., if 



Vp = + ^ < & ' « yV? + ^ + (1/12)0 /2 + (l/12)A 2 e- 2 ^- 6fi . (5.12) 

30 



Clearly, this condition is fulfilled in every case unless 0' = A = (no dilaton and axion - that 
is, the general relativity vacuum regime), which reflects the fact that a particle cannot be 
scattered infinitely many times and that there is no chaos in the Einstein frame. This result 
is expected since in the Einstein frame both dilaton and axion fields behave as stiff fluids 
with the equation of state p = g. A numerical dicussion of 10- dimensional axion-dilaton 
low-energy effective-action models in the Einstein frame where nine dimensions were split 
into three isotropic 3-dimensional spaces leading effectively to our anisotropic 4-dimensional 



model, was also given in Ref. |25| with the same final conclusion about the non-existence of 
chaos. 

B. String Frame and Axion frame 

In the string frame we can use the same parametrization as in (5.1), but we just drop the 
tildes. The potential (5.2) can also be used without tildes. In that parametrization Eq.(3.5) 
becomes 



a 2 = ii + # - \ft + \e~ 2a V^ ± ) + ±-A 2 e- ea + aj>. (5.13) 
oo 12 



After applying the variables (3 and <ft defined by (4.5) we can remove the <fra term, obtaining 
f = f + Q^p 2 + + 6^! + e~ 2 ^V(ip±) + ^A 2 e~ 2 ^. (5.14) 



Following the analysis given in ||15|| , we apply a new time coordinate r defined as 

dt = dre-*, (5.15) 

and define a new variable, y, which is the logarithm of an averaged scale factor in the 
conformally related axion frame (or B-frame) in which the axion is minimally coupled |24 . 
This is given by 

y = V3<f) + (3 — —= In abc, (5. 16) 

v3 

and brings (5.14) to the form ((...) T = d{...)/dr) 
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y 2 = <j? T + 12^ r + u^_ t + A 2 e- 2<f > + 2e-7SV(V>±). 



(5.17) 



Eq. (5.17) is, in fact, the Hamiltonian constraint obtained from the action (A^ 



^ttG) 
(5.18) 



The canonical momenta are then 



7I> 



A, 



7T+ = 6A s ^+r, 
7T_ = 6A s ^-r, 



and the Hamiltonian is just 



T " 2 ~ n v ' 12 



(5.19) 



(5.20) 



4 2 

Following the analysis of the previous subsection V.A., we see that the maximum apparent 
velocity of the wall at ip + <C — 1, ty- ~ is given by 

1 



Vr 



max 2V3 

and the condition for chaotic scatterings to cease is just that 



(5.21) 



(5.22) 



12 TT 12 

This is clearly fulfilled except in the general relativity case where = A = (i.e., no axion 
and dilaton fields). This gives our final conclusion that there is no chaos in BIX string 
cosmology in the string or axion frames. 

It is not surprising that the physical behaviour should be similar in every frame ||, so 
that if there is no chaos in the Einstein frame there should not be chaos in any other frame. 
String theory appears to impose too much symmetry through its duality invariances for 
chaos to appear Pi 



10 Chaos has been studied in other related situations in Ref. |26| 
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VI. DISCUSSION 



In this paper we have carried out a detailed analysis of the spatially homogeneous uni- 
verses of Bianchi type IX in the context of low-energy string theory. These universes are of 
special interest in relativistic cosmology because they display chaotic behaviour in vacuum 
and in the presence of fluids with p < p. They were originally termed 'Mixmaster' universes 
by Misner because they offered the possibility for light to travel all the way around the 
universe in different directions. Moreover, they are the most general closed universes which 
are spatially homogeneous and may be closely related to parts of the general solution of 
Einstein's equations in the neighbourhood of a strong curvature singularity of the sort that 
characterises the initial state of general relativistic cosmological models. This behaviour has 
also been extensively investigated because it is of intrinsic mathematical interest. It is also 
known that its occurrence in general relativity depends upon the dimensionality of space. 
We investigated the string cosmological equations for the type IX metric. We found that 
chaotic behaviour does not occur in string cosmology in either the Einstein or the string or 
axion frames. While it is possible for finite sequences of oscillations to occur in the scale 
factors' evolution on approach to t — 0, these oscillations cannot continue indefinitely. They 
inevitably terminate in a state in which all the three orthogonal scale factors decrease with 
decreasing time monotonically on approach to the initial singularity. We investigated the 
detailed sequences of evolutionary changes that can take place in the evolution during the 
finite sequences of oscillations between epochs which are well approximated by Kasner uni- 
verses. We found that the duality symmetry required of the string evolution introduced new 
invariances for the possible changes in the Kasner parameters in addition to those which 
characterise the Kasner-to-Kasner cycles of oscillations. The requirements of duality invari- 
ance on the evolution of the metric appear to be so constraining that chaotic behaviour is 
excluded. We have obtained these results in two complementary ways: by direct matching 
of asymptotic expansions to the solutions of the sytem of non-linear ordinary differential 
equations of string cosmology and by use of the Hamiltonian formulation of cosmology. In 
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the Hamiltonian picture the evolution of the type IX string cosmology is represented as the 
motion of a 'universe point' inside a potential that is open only along three narrow chan- 
nels. The walls of this 'almost closed' potential of Henon-Heiles type expand outwards as 
the singularity is reached. Whereas in the vacuum models of general relativity, the universe 
point always catches the walls and bounces chaotically around within the potential, in string 
theory the universe point need never catch the walls. If it is moving towards a wall at a 
very oblique angle then the normal component of its velocity towards the wall can become 
too small for it ever to catch the wall. In general, we find that this situation always arises 
after a finite number of collisions have occurred in the Mixmaster string cosmology. The 
resulting asymptotic state is therefore similar to that in a model with no potential walls at 
all; that is the Bianchi type I or Kasner universe. 
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APPENDIX A: RICCI TENSOR IN A COORDINATE FRAME 

In this Appendix we discuss a possibility of admittinig a time-dependent antisymmetric 
tensor potential B^ v = B^it) as given in Ref. || to the axisymmetric Bianchi IX model. 
It has been proven [|J that such a potential cannot be admitted to a general (i.e. non- 
axisymmetric case). However, to achieve this, we first give the components of Ricci tensor 
for the axisymmetric Bianchi IX model in terms of coordinates rather than in orthonormal 
frames of Section II. For the sake of generality, we start with a general metric. 

The metric (2.9) in a coordinate frame has the following components 
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9m = i, 

911 = ~c 2 (t), 

a 2 (t) cos 2 ip + b 2 (t) sin 2 ip 



922 
933 



sin 2 9 



a 2 {t) sin 2 if) + b 2 (t) cos 2 if) 
913 = -c 2 (t) cos 0, 

923 — — sin -0 cos if> sin 9 a 2 (t) — b 2 (t) 



c 2 {t) cos 2 9, 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 
(A6) 



The Ricci tensor components in the coordinate frame can be calculated using the relations 
H (cf. Eqs. (2.9)-(2.12)) 



Ra/3 — ^a^pR-jki 
Tlf3 _ j 3 r>k 



(A7) 
(A8) 



where R a p is the Ricci tensor in the coordinate frame while Rjk is the Ricci tensor in the 
orthonormal frame (correspondingly a, (3 = 0,1,2,3 are the coordinate frame indices and 
i,j = 0,1,2,3 are the orthonormal frame indices). Thus, 1,2,3 refer to a 1 , a 2 , a 3 , while 
1, 2, 3 refer to ip, 9, respectively. According to (2.10)-(2.12), 

e\ = a cos ip, 
e\ = a sin ip sin 9, 
e\ = bsintp, 
e| = — bcosijj sin 9, 
ef = c, 
e| = ccos^, 



(A9) 



and 



sin if) cot 9, 



(A10) 



- cos if) cot 9, 
b 
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3 ~ c' 



e 1 — — cos ■?/>, 
a 

e 2 = -sm^, 

3 1 sin if) 
1 a sin # ' 

3 1 COS if) 



(All) 



2 6sin#' 

Then, the Ricci tensor components in the coordinate frame are given by 
Rii = c 2 R 33 , 

i?25 — ° 2 cos2 ^-Rn + & 2 sin 2 1PR22, 

i?33 = a 2 sin 2 if) sin 2 6*-Rn + b 2 cos 2 ^ sin 2 6*-R 2 2 + c 2 cos 2 6^33, 
i?i3 = c 2 cos9R 33 , 

R23 = sin t/> cos ^ sin 9 {a 2 Rn — b 2 R 2 2 
0. 



and 



R 



12 



4 



4 

Rl 



cos 2 ipR{ + sin 2 V'-R 2 , 
sin 2 if)R\ + cos 2 V'-R-L 
- sin 2 if)R\ - cos 2 + #3 
sin if) cos if) sin # (i?} — R\ 



cos 6> 



i?2 



1 . 



i?2 
R\ 



sin ^ cos ^ 
sin 9 

— sin ^ cos if) cot 6* (i?J — R\ 

0, 
0. 



(A12) 
(A13) 
(A14) 
(A15) 
(A16) 
(A17) 

(A18) 
(A19) 
(A20) 
(A21) 
(A22) 
(A23) 
(A24) 
(A25) 
(A26) 



If two axes are the same {a{t) = b(t)) the metric (2.9) (or its components given by 
(A1)-(A6)) simplifies to 0,[ll 
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ds 2 = dt 2 - c 2 + cos6dip) 2 - a 2 (d9 2 + sin 2 9d<p 2 ) . 
The nonvanishing Christoffel symbols for the metric (A27) are 



(A27) 



1 03 



r 1 --- 

01 ~ c ' 

c d\ 

cos fl, 

c a/ 



1 02 — 



r?i = cc, 



1 03 — 



J. 22 — aa 



r§3 = cc cos 2 # + da sin 2 9, T j 5 = cc cos 6* 



1 21 



lc 2 
2 a 2 



cot 0, r| 



1 



23 



2sin# 



cos 2 9 - 1 



lc 2 
2 a 2 



, - —^sm.0, r§g = sin 9 cos 9 



c 2 -a 2 



(A28) 



V-- — — - 

1 12 _ 



lc 2 



/ 1 c 2 

= cot^ l--- 

\ 2 a" 1 



2 a 2 sin 

The gradients of the dilaton calculated with respect to the metric (A27) are given by 

(A29) 
(A30) 
(A31) 
(A32) 



V 5 V°0 = 


0, 


VivV = 


C 


V 5 V 2 = 


V S V 5 = 

a 


VgVV = 


<j) ( - — — ) cos#, 
Vc a/ 



and 



V o 0VV = 2 



(A33) 



The nonzero components of the Ricci tensor in a coordinate frame for the metric (A27) 
can be obtained from (A12)-(A26) by putting R\ = R\ and using (2.17)-(2.20) for a = b or 
directly using the Christoffel symbols. This gives 



d lc 4 

Rll = CC + 2(X h -— r, -R13 = i?n cos#, 

a 2 a 4 

c 1 c 2 C a 

R22 = da + d 2 + da- + 1 -, _Rgo = 2-, 

c 2 a 2 c a 

i?33 = R n cos 2 + R22 sin 2 0, 



(A34) 



so 
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-^ = - + 2-, (A35) 
c a 

= -+2-- + £4. (A36) 
c a c 2 a 4 

o 5 a d 2 dc 1 1 c 2 . , 

-i^ = -i^ = - + — + — + — - --, A37 
d a a 2 a c a 2 2 a 4 

n i (a c ac a 2 1 c 2 \ . ,. no , 

+ — + — -— cosfl, (A38) 

\a c ac a z a z or J 

and the Ricci scalar reads 

c a a c a 2 1 1 c 2 . . 

i2 = _ 2 - - 4- - 4— - 2- - 2- + -- (A39) 
c a ac cr a z 2 a 4 

Note that in the coordinate frame R 2 . = while in the orthonormal frame (cf. Eq.(2.18)- 
(2.19)) R\ = R\. This is reasonable, since the metric tensor is given by (A27), for the former 
case, and by (2.9) with a = b for the latter case, where the indices refer to the orthonormal 
basis rather than to the chosen coordinates. 

Now we convert our notation in terms of the three-index torsion field H to the notation 
given in using the pseudoscalar torsion field h. Following || we define (a, (3-,n,v = 
0,1,2,2) 

= eV^V. (A40) 

The equation of motion for /i-field obtained via integrability conditions is then 

VV^/i + V M 0V M /i = 0. (A41) 

One can easily see that for the time-independent antisymmetric tensor potential B^ v = 
B^ u (x) the h field can only depend on time and the equation (A41) reads as 

h+(- + l + -)h + j>h = 0, (A42) 



which integrates to give 



so from (A40) we have 



a b c , 



h = -A*——, (A43) 
abc 
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or 



#123 = AsinO, (A45) 

as required by (2.25) and H 2 = —QA 2 /a 2 b 2 c 2 . With the H field chosen as above the equation 
of motion (2.3) is easily fulfilled. There is also a trivial solution of (A42), h = 0, but it 
corresponds to a constant torsion field. For the time-independent pseudoscalar axion field 
h = h(x) (time-dependent antisymmetric tensor potential = B^ u (t)) the equation of 
motion (A41) reads 

d^d^h + T^d p h = 0, (A46) 

which for the metric (2.9) reads as 

g ll d\h + g 22 d\h + g m d\h + g^d-^h + g 22 cot 9d 2 h = 0. (A47) 

For simplicity, let us introduce 

(H;f = H^H"* = -2e 2 * (fyf - Stf") d £ hd P h, (A48) 
H 2 = H mP H^. (A49) 

The non-zero components of these quantities (the energy- momentum tensor) which are used 
in the field equations (2.1)-(2.3) for the metric components (A27) are given by 

(jf°) 2 = -2e 2 ^ (g ll dihdih + 2g r3 dihd 3 h + g 22 d 2 hd 2 h + g m d 3 hd 3 h) , (A50) 

{H\f = -2e 2 * (g l %kdsh + g 22 d 2 hd 2 h + g m d 3 hd 3 h) , (A51) 

(H~l) 2 = -2e 2(j} (g ll dihdih + 2g r3 dihd 3 h + g m d 3 hd 3 h) , (A52) 

(#|) 2 = -2e 2 ^ (g ll dihdih + g r3 dihd 3 h + g 22 d 2 hd 2 h) , (A53) 

(Hi) 2 = 2e 2 * (g^hd-sh + g l3 d 3 hd 3 h) , (A54) 



and 
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H 2 = -6e 2 ^ (g ll dihdih + 2g I ~ 3 d- 1 hd 3 -h + g 22 d- 2 hd- 2 h + g 3 ~ 3 d- 3 hd- 3 h) , (A55) 
With the choice 

d 3 h = E sin 6 ^ ^/i = = 0, (A56) 

the field equations remain homogeneous and 



(^) 2 = (H{? = {Hlf = 2^-, (A57) 



E 2 e 2* 

(Hi) 2 = 0, (A58) 

E 2 e 2 <t> 

(Hi) 2 = 2=4- cos ^, (A59) 

H 2 = 6^— (A60) 
a 2 

and the equation of motion (A47) becomes 

g™d v (Enn9) = 0, (A61) 
and it is satisfied. Finally, from (A40), we have 



H 012 = (A62) 



Ee* 
crc 

H m = Ece+, (A63) 



and H 2 = 6Ee 2< ^/a 2 . This obeys the axion equation of motion (2.3) but is in contradiction 
with the axisymmetry condition for the Ricci components (A37) since R 2 = R 3 there while 
(Hi) 2 = and (iff) 2 ^ here. 

A possible ansatz which would fulfil the axisymmetry condition would be 

Ba(t)sm9 IK . 

dih= . y (A64) 

^/a(*) 2 sin 2 + c(i) 2 cos 2 
but it leads to both time and space dependences of the pseudoscalar axion field, h = 
h(t,ip,6), and does not fulfil the equation of motion (A41) (nor (A47) which is obtained 
for h = h(x)). 

One could also try to add 
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d- 2 h = D = const., (A65) 

to the nonzero component (A56), but again this does not fulfil the equation (A47). 

The final conclusion is that one is not able to impose the axion field even in axisymmetric 
BIX models despite the fact that there is a distinguished direction in the model (which is 
different from electromagnetic field case - see [IT]). The reason seems to be that even in the 



axisymmetric case there is still SO(3) symmetry group present and we are only adding an 
additional symmetry SO(2) which does not cancel the former one, giving the total symmetry 
SO(3) (8)50(2) rather than just S0{2). 
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FIGURES 

FIG. 1. The regions I, II, III of permissible values of the parameters u and M for the 
parametrization (3.60) of the Kasner indices p\,P2 and p%. We restricted ourselves to — 10 < u < 10 
although regions I and III extend to infinity. Special isotropic FRW (+) and FRW (-) points are 
given for (a/3 - a/2, -1 + a/3) and (-a/3 + a/2, -1 - a/3) respectively. 

FIG. 2. The regions I, II, III of permissible values of the parameters u and M for the 
parametrization (3.62) of the Kasner indices p\ , p2 and p%. We restricted ourselves to — 10 < u < 10 
although regions I and III extend to infinity. Special isotropic FRW (+) and FRW (-) points are 
given for (-a/3 + a/2, -1 + v^) and (a/3 - a/2, -1 - a/3) respectively. 

FIG. 3. The plot of the conditions (3.96)-(3.97) for the transitions from one Kasner epoch 
to another to begin in terms of the first and the "fourth" Kasner indices p\ and M. Obviously, 
isotropic FRW (+) and FRW (-) points (-1 + a/3, 1/a/3) and (-1 - a/3, -1/a/3) respectively, are 
excluded. The transitions occur below the line p = M/2. From the picture one can easily see that in 
vacuum general relativity case (M = 0) the transitions are possible whenever — 1/3 < p\ < 0. This 
is in agreement with the standard calculations. However, in the dual case (M = —2) transitions 
do not occur at all. 

FIG. 4. The plot of the conditions (3.101)-(3.103) for the transitions from one Kasner epoch 
to another to begin, in terms of the two Kasner indices p\ and p2- As in Fig. 3 the isotropic FRW 
(+) and FRW (-) points, (1/a/3, 1/a/3) and (— 1/a/3, — 1/a/3) respectively, together with the two 
neighbouring regions, are excluded. 
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Figure 1: 
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Figure 2: 
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Figure 3: 




48 



Figure 4: 
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